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Abstract 

> 

We calculate the nucleon sigma term in two-flavor lattice QCD utilizing the Feynman-Hellman 

^t- ' 

theorem. Both sea and valence quarks are described by the overlap fermion formulation, which 
\Q . preserves exact chiral and flavor symmetries on the lattice. We analyse the lattice data for the 

o , 

nucleon mass using the analytical formulae derived from the baryon chiral perturbation theory. 



From the data at valence quark mass set different from sea quark mass, we may extract the 



^ ■ sea quark contribution to the sigma term, which corresponds to the strange quark content. We 

H \ 

find that the strange quark content is much smaller than the previous lattice calculations and 
phenomenological estimates. 

PACS numbers: 
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I. INTRODUCTION 



A piece of information on the nucleon structure can be extracted from its quark mass 
dependence. Nucleon sigma term a^N characterizes the effect of finite quark mass on the 
nucleon mass. Up to non-analytic and higher order terms, the nucleon mass is written as 
M/v = Mq + (Tttn, where Mq is the nucleon mass in the chiral limit. The exact definition of 
ovtv is given by the form of a scalar form factor of the nucleon at zero recoil as 

ct-kN = m ud ((N\uu + dd\N) - V(0\uu + dd\0)) (1) 

where m u d denotes degenerate up and down quark mass. The second term in the parenthe- 
sis represents a subtraction of the vacuum contribution, and V is the (three-dimensional) 
physical volume 1 . For the sake of simplicity we represent the vacuum subtracted matrix 
element (N\qq\N) — V(0\qq\0) by (N\qq\N) in what follows, (q represents a quark field: up 
(u), down (d), or strange (s).) Note that the sigma term is renormalization group invariant, 
since the renormalization factor cancels between the quark mass m q and the scalar operator 
qq. 

While the up and down quarks contribute to a^N both as valence and sea quarks, the 
strange quark appears only quark contribution. As a measure of the strange quark 

content of the nucleon, the y parameter 

v = 2(N\ss\N) 

y {N\uu + dd\NY y ' 

is commonly introduced. Besides characterizing the purely sea quark content of the nucleon, 
which implies a clear distinction from the quark model picture of hadrons, this parameter 
plays an important role to determine the detection rate ofpossible neutralino dark matter 
in the supersymmetric extension of the Standard Model 

aaaaaaa 

. Already with the 

present direct dark matter search experiments one may probe a part of the MSSM model 
parameter space, and new experiments such as XMASS and SuperCDM will be able to 
improve the sensitivity by 2-3 orders of magnitude. Therefore, a precise calculation of the 
y parameter (or equivalently another parameter fx s = m s {N\ss\ N)/M^) will be important 
for excluding or proving the neutralino dark matter scenario. 



1 The nucleon state \N{p)) is normalized as {N(p)\N{p')) — (27r) 3 S^ (jp— p'). In (TjJ we omit the momentum 
argument for the nucleon, since we do not consider finite momentum insertion in this paper. 
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Phenomenologically, the sigma term can be related to the ttN scattering amplitude at a 



certain kinematical point, i.e. the so-cal 
is in the range T, C d = 70 ~ 90 MeV 



ed Cheng-Dashen point t = +2m^ [8|. Its value 
j. After the corrections for the finite value of t, 
which amounts to —15 MeV [10], one obtains a n N = 55 ~ 75 MeV. On the other hand, 
the octet breaking of the nucleon mass, or the matrix element (N\uu + dd — 2ss\N), can 
be evaluated from the baryon mass spectrum. At the leading order of Chiral Perturbation 
Theory (ChPT), the value of the corresponding sigma term is a ~ 26 MeV, while the heavy 
baryon ChPT (BChPT) gives a = 36 ± 7 MeV ll|. The difference between g^n an d o is 
understood as the strange quark contribution; algebraically the relation is a n N = c*7(l — y)- 
Then, one obtains a large value of y: y=0.3-0.6. (The value of y is even larger than the 



estimate y ~ 0.2 in 



lOj, because of the more recent experimental data [9].) For other 



phenomenological estimates, see, e.g. [12j |. Such large values of y cannot be understood 
within the valence quark picture, hence raises a serious problem in the understanding of the 
nucleon structure. We note however that the analysis within chiral effective theories suffers 
from significant uncertainties of the low energy constants, especially at higher orders. 

Using lattice QCD, one can in principle calculate the nucleon sigma term without in- 
volving any model parameters, since lattice calculation for a wide range of quark masses in 
the chiral regime offers essential information on the low energy constants which cannot be 
determined by experimental data alone. Furthermore, it is possible to determine the valence 
and sea quark contributions separately. A direct method to extract them is to calculate 
three-point functions on the lattice including an insertion of the scalar operator. It can also 
be done in an indirect way by analyzing the quark mass dependence of the nucleon mass 
for valence and sea quarks separately. Obviously, the dynamical fermion simulations are 
necessary to extract the disconnected contributions in the indirect method. 

Previous lattice results were g^n = 40-60 MeV, y = 0.66(15) [IJbJ], and a n N = 50(3) MeV, 
y = 0.36(3) jl^| within the quenched approximation, while a two-flavor QCD calculation [ljj] 
gave a n N = 18(5) MeV and y = 0.59(13). There are apparent puzzles in these results: firstly 
the strange quark content due to the disconnected diagram (the value of y) is unnaturally 
large compared to the up and down contributions that contain the connected diagrams too. 
Secondly the values of the sigma term in the quenched and unquenched calculations are 
rather different, which might also imply significant effects of quark loops in the sea. 

Concerning the first point, it was pointed out that using the Wilson-type fermions, which 
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violate the chiral symmetry on the lattice, the sea quark mass dependence of the additive 
mass renormalization and lattice spacing can give rise to a significant uncertainty in the 
sea quark content 16] . Unfortunately, after subtracting this contamination the unquenched 
result has large statistical error, y = —0.28(33). In the present work, we remove this 
problem by explicitly maintaining exact chiral symmetry on the lattice for both sea and 
valence sectors, as described below. 

The second puzzle may be resolved by incorporating an enhancement due to pion loops. 
Within BChPT at 0(p 3 ) or 0(p 4 ), a curvature is expected in the quark mass dependence 
of the nucleon, hence the sigma term, a derivative of in terms of m q , increases towards 
the chiral limit. 

An analysis using existing lattice data of two-flavor QCD with m n > 550 MeV by CP- 
PACS 12|, JLQCD 18], and QCDSF 19] yields a wN = 48 ± MeV Q, 21], which is 



slightly smaller than but is still consistent with the phenomenological analysis. A more 
recent lattice data by the ETM Collaboration with m n = 300-500 MeV in two-flavor QCD 
reported a higher value a^N = 67(8) MeV [22]. Such an analysis for the disconnected 
contribution to extract the strange content is yet to be done, which is another main point 
of this work. 

In this work, we analyze the data of the nucleon mass obtained from a two-flavor QCD 



simulation employing the overlap fermion 



231 ] . (For other physics results from this simu- 



lation, we refer 24j and references therein.) The overlap fermion 25j, |26J preserves exact 
chiral symmetry on the lattice, and there is no problem of the additive mass shift of the 
scalar density operator, that was a main source of the large systematic error in the previous 
calculations of the sigma term. We use the overlap fermion to describe both the sea and 
valence quarks. Statistically independent ensembles of gauge configurations are generated 
at six different sea quark masses; the nucleon mass is measured for various valence quark 
masses on each of those gauge ensembles. Therefore, we are able to analyze the valence and 
sea quark mass dependence independently to extract the connected and disconnected contri- 
butions. An estimate of the strange quark content can thus be obtained in two-flavor QCD. 
In the analysis, we use the partially quenched BChPT, which corresponds to the lattice cal- 
culations with valence quark masses taken differently from the sea quark masses. Therefore, 
the enhancement of the sigma term towards the chiral limit is incorporated. Since the two- 
flavor QCD calculation cannot avoid the systematic error due to the neglected strange sea 
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quarks, our result should not be taken as a final result from lattice QCD. Nevertheless our 
study with exact chiral symmetry reveals the underlying systematic effects in the calculation 
of the nucleon sigma term, especially in the extraction of its disconnected contribution. It 
therefore provides a realistic test case, which will be followed by the 2+1-flavor calculations 
in the near future 2 . 

Our paper is organized as follows. In Section [Til we introduce the basic methods to 
calculate the nucleon sigma term. Our simulation set-up is described in Section IIII1 Then, 
in Section HVl we describe the BChPT fit to obtain the sigma term. In Section [V] we study 
the sea quark content of the nucleon from PQChPT. In Section Wl\ we compare our results 
with previous calculations and discuss the origin of the discrepancy. Our conclusion is given 
in Section IVTT1 



II. METHOD FOR CALCULATING NUCLEON SIGMA TERM 

The matrix element defining the nucleon sigma term (CD) can be related to the quark mass 
dependence of the nucleon mass using the Feynman-Hellman theorem. Consider a two-point 
function of the nucleon interpolating operator Ojv(£, x) 

G(t) = f d 3 x(0\O N (t,x)O* N (0,0)\a) = Z- 1 j VA^ l[(VqVq)O N (t,x)Ol(0,6)e- s (3) 

J J q 

with the QCD action S defined by the gluon field strength F„ v and the quark field q as 



(4) 



and the partition function Z. The sum in (j3J) runs over flavors (q = u and d) according to 
the underlying two-flavor theory (Nf = 2). By taking a partial derivative with respect to a 
valence quark mass m val or a sea quark mass m sea corresponding to the degenerate u and d 



For a very recent result from 2+1-flavor QCD, see 27|. 
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quark masses m u d (= m u = rrid), we obtain 



dGjt) 
dm val 

dGjt) 
<9m SPa , 



- / d 3 x ( 



- / d 3 x ( 



+G(t) ( 



0^,^)0^(0,0) 
0^,^)0^(0,0) 



d 4 y (^)(y) 

q=u,d 



d A y (m)(v) 



q=u,d 



(5) 



d 4 y Y M (y) 

q=u,d 



disc 



(6) 



The subscripts "conn" and "disc" on the expectation values indicate that only the connected 
or disconnected quark line contractions are evaluated, respectively. 

Dividing the integration region of t y , a temporal component oiy, into three parts, i. e. t y < 
0, < t y < t, and t < t y , and inserting the complete set of states between the operators, one 
can express G(t), dG(t)/dm va \, and dG(t)/dm sea in terms of matrix elements. Comparing 
the leading contribution at large t behaving as texp (—M^t) with the nucleon mass, we 
obtain the relations 



dM, 



N 



dm va i 
dM N 



(N\(uu + dd)\N) conn 
(N\(uu + dd)\N) disc , 



(7) 
(8) 



Note that the short-hand notation to omit the term — V (0\(qq)\0) applies only for the dis- 
connected piece. 

This derivation of the Feynman-Hellman theorem does not assume anything about the 
renormalization scheme nor the regularization scheme. The contact terms in (jSJ) and (jUJ) 
are irrelevant for the formulas ([7D and (|5J), since only the long-distance behavior of the 
correlators is used. 

In the present study we exploit this indirect method to extract the matrix elements 
corresponding to the nucleon sigma term. 

Another possible method to calculate the nucleon sigma term is to directly calculate the 
matrix element from three-point functions with an insertion of the scalar density operator 
(uu + dd)(x), as carried out, e.g. in Q, Q in the quenched approximation. In principle, it 
gives a mathematically equivalent quantity to the indirect method, provided that the indirect 
method is applied with data at sufficiently many sets of {m va [,m sea ) so that the derivatives 
are reliably extracted. The order of the derivative and the path integral does not make any 



difference at finite lattice spacing and volume. Numerical difference could arise only from 
the statistical error and the systematic error in the fit of the data. A practical advantage of 
the indirect method is that the sum over the position of (uu + dd) (x) is automatic, whereas 
in the direct method it must be taken explicitly to improve statistical accuracy. On the 
other hand, the direct method is more flexible, as one can take arbitrary quark masses for 
the "probe quark" to make a disconnected loop from the (uu + dd)(x) operator, while in 
the indirect method the probe quark mass is tied to the sea quark mass. Therefore, we can 
only estimate the strange quark content from the calculation done at the strange quark mass 
equal to the sea quark mass as dM N /dm B ^\ mval=nhm = ma = 2(iV|ss|iV} disc . 



III. LATTICE SIMULATION 



We make an analysis of the nucleon mass using the lattice data obtained on two-flavor 
QCD configurations generated with dynamical overlap fermions 23J]. The lattice size is 
16 3 x 32, which roughly corresponds to the physical volume (1.9 fm) 3 x(3.8 fm) with the 
lattice spacing determined through the Sommer scale rn as described below. The overlap 



fermion is defined with the overlap-Dirac operator 



25 



26] 



/ Tfl \ f Tfl \ 

D(m q ) = [m + -^J + (m - -yj 75Sgn [H w (-m )} (9) 

for a finite (bare) quark mass m q . The kernel operator Hw(-mo) = 7s£V( — m ) is con- 
structed from the conventional Wilson-Dirac operator Dy/(— m ) at a large negative mass 
—niQ. (We set mo = 1.6 in this work.) For the gluon part, the Iwasaki action is used at (3 
= 2.30 together with unphysical Wilson fermions and associated twisted-mass ghosts 28] 
introduced to suppress unphysical near-zero modes of Hw(— With these extra terms, 
the numerical operation for applying the overlap-Dirac operator (Q is substantially reduced. 
Furthermore, since the exact zero eigenvalue is forbidden, the global topological charge Q 
is preserved during the molecular dynamics evolution of the gauge field. Our main runs are 
performed at the trivial topological sector Q — 0. For each sea quark mass listed below, 
we accumulate 10,000 trajectories; the calculation of the nucleon mass is done at every 20 
trajectories, thus we have 500 samples for each m sea . For more details of the configuration 
generation, we refer to 231 ] . 

For the sea quark mass am sea we take six values: 0.015, 0.025, 0.035, 0.050, 0.070, and 
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am sea a [fm] m w [GeV] 


m^L 


0.015 0.1194(15) 0.2880(18) 


2.8 


0.025 0.1206(18) 0.3671(13) 


3.5 


0.035 0.1215(15) 0.4358(13) 


4.2 


0.050 0.1236(14) 0.5217(13) 


5.0 


0.070 0.1251(13) 0.6214(11) 


6.0 


0.100 0.1272(12) 0.7516(14) 


7.2 



TABLE I: Lattice spacing and pion mass calculated for each sea quark mass. 



0.100 that cover the mass range m s /6-m s with m s the physical strange quark mass. Analysis 
of the pion mass and decay constant on this data set is found in [29J. 

The lattice spacing determined through the Sommer scale ro of the static quark poten- 
tial slightly depends on the sea quark mass; the numerical results are listed in Table [B 
Extrapolating to the chiral limit, we obtain a = 0.118(2) fm assuming the physical value 
r = 0.49 fm. In the following analysis, we use this value to convert the lattice results to 
the physical unit. 

The two-point functions, from which the nucleon mass is extracted, are constructed from 
quark propagators described by the overlap fermion. In order to improve the statistical 
accuracy, we use the low-mode preconditioning technique [30j], i.e. the piece of the two-point 
function made of the low modes of the overlap-Dirac operator is averaged over many source 
points. In our case, the source points are set at the origin on each time slice and averaged 
over different time slices with 50 chiral pairs of low modes. For the source to solve the quark 
propagator, we take a smeared source defined by a function 0(|x|) oc exp(— A\x\) with a 
fixed A = 0.40. We then calculate the smeared- local two-point correlator and fit the data 
with a single exponential function after averaging over forward and backward propagating 
states in time. The statistical error is estimated using the standard jackknife method with 
a bin size of 10 samples, which corresponds to 200 trajectories. In the calculation of the 
nucleon mass, we take the valence quark masses am va i = 0.015, 0.025, 0.035, 0.050, 0.060, 
0.070, 0.080, 0.090, and 0.100. 

Figure [1] shows an effective mass plot of the nucleon at am sea = 0.035. Data are shown for 
the nucleon made of degenerate valence quarks at the nine available masses. We find a good 
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FIG. 1: Effective mass of the smeared-local nucleon correlator. Data are shown for various degen- 
erate valence quark masses m va i at a fixed sea quark mass am sea = 0.035. 

plateau for all sea and valence quark mass combinations; the fit with a single exponential 
function is made in the range [5,10]. The fitted results are shown by thick horizontal lines 
in Figure [1] and summarized in Table [III 

IV. ANALYSIS OF THE UNITARY POINTS WITH BARYON CHIRAL PER- 
TURBATION THEORY 

A. Naive fits with BChPT 



In this section, we analyze the lattice data taken at the unitary points, i. e. sea and valence 
quarks are degenerate. In this case, the conventional baryon chiral perturbation theory 
(BChPT) 3l| for two flavors is a valid framework to describe the quark mass dependence of 
the nucleon. It develops a non-analytic quark mass dependence and leads to the enhancement 
of the nucleon sigma term near the chiral limit. 

In BChPT, the nucleon mass is expanded in terms of the light quark mass or equivalently 
pion mass squared m^. We follow the analysis done in 20]. The expression for the nucleon 
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am sca am val 


am ps 


aM N 


am sca am val 


amps 


aM N 


0.015 0.015 


0.1729(12) 


0.6647(59) 


0.050 0.015 


0.17402(87) 


0.6895(58) 


0.025 


0.2210(10) 


0.7038(47) 


0.025 


0.22279(82) 


0.7274(43) 


0.035 


0.25998(95) 


0.7381(44) 


0.035 


0.26218(80) 


0.7606(39) 


0.050 


0.30966(93) 


0.7858(43) 


0.050 


0.31228(80) 


0.8072(38) 


0.060 


0.33918(96) 


0.8164(43) 


0.060 


0.34199(82) 


0.8369(40) 


0.070 


0.3668(10) 


0.8460(44) 


0.070 


0.36967(84) 


0.8657(39) 


0.080 


0.3929(11) 


0.8750(45) 


0.080 


0.39583(93) 


0.8937(41) 


0.090 


0.4180(12) 


0.9033(47) 


0.090 


0.4208(10) 


0.9211(43) 


0.100 


0.4423(15) 


0.9312(48) 


0.100 


0.4449(13) 


0.9480(45) 


0.025 0.015 


0.17185(99) 


0.6597(60) 


0.070 0.015 


0.17512(76) 


0.6870(63) 


0.025 


0.21990(81) 


0.6960(45) 


0.025 


0.22444(65) 


0.7259(44) 


0.035 


0.25906(76) 


0.7300(40) 


0.035 


0.26399(62) 


0.7610(39) 


0.050 


0.30900(72) 


0.7780(37) 


0.050 


0.31414(64) 


0.8098(37) 


0.060 


0.33864(72) 


0.8084(37) 


0.060 


0.34388(68) 


0.8407(37) 


0.070 


0.36625(75) 


0.8380(39) 


0.070 


0.37166(75) 


0.8705(38) 


0.080 


0.39230(81) 


0.8667(38) 


0.080 


0.39800(88) 


0.8996(39) 


0.090 


0.41711(92) 


0.8946(39) 


0.090 


0.4232(11) 


0.9280(41) 


0.100 


0.4409(11) 


0.9220(42) 


0.100 


0.4477(14) 


0.9559(42) 


0.035 0.015 


0.17299(96) 


0.6859(73) 


0.100 0.015 


0.17663(71) 


0.7040(65) 


0.025 


0.22168(84) 


0.7186(45) 


0.025 


0.22563(61) 


0.7419(44) 


0.035 


0.26111(81) 


0.7505(35) 


0.035 


0.26548(58) 


0.7761(37) 


0.050 


0.31136(82) 


0.7966(30) 


0.050 


0.31619(56) 


0.8236(35) 


0.060 


0.34124(85) 


0.8263(30) 


0.060 


0.34621(57) 


0.8536(34) 


0.070 


0.36917(90) 


0.8553(30) 


0.070 


0.37418(61) 


0.8827(34) 


0.080 


0.39567(99) 


0.8837(30) 


0.080 


0.40061(70) 


0.9110(35) 


0.090 


0.4211(11) 


0.9115(32) 


0.090 


0.42587(87) 


0.9387(36) 


0.100 


0.4456(14) 


0.9388(33) 


0.100 


0.4502(11) 


0.9659(37) 



TABLE II: Numerical results for the pseudo-scalar meson mass and the nucleon mass for each sea 
and valence quark masses. 
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mass to the order 0(p 3 ) has a form 
M N = M - 4 Cl m 2 - -^^ml 4 



e£(/x) - 



3<? 



647r 2 ^Mo 



1 + 2 log 



256vr/2M 2 
and that to 0(p 4 ) is 



(10) 



M, 



N 



Mn 



+ 



4cim 2 



3^ 2 



327T/2 



1 mi 



3^ 



c 2 



647r 2 /|V^"o 2; 32vrV 2 V^o 



9a 



- 8ci + c 2 + 4c; 



log — 



256tt/2M ( ; 



(11) 



There are many parameters involved in these expressions. First of all, M is the nucleon 
mass in the chiral limit and f n is the pion decay constant. The constant qa describes the 
nucleon axial-vector coupling. Its experimental value determined by the neutron (3 decay 
is 9a — 1-270(3) 32j. The parameters ci, c 2 , and C3 are low energy constants (LECs) at 
0(p 2 ); their phenomenological values are c\ = — 0.9~^q 2 GeV -1 , c 2 = 3.3 ± 0.2 GeV -1 , and 
C3 = —A.7t\'x> GeV -1 (for a summary, see, for example [3]). In the fit using (TTTi) dis- 
cussed below, we fix (c 2 ,c 3 ) at two representative combinations, (3.2 GeV -1 ,— 3.4 GeV -1 ) 
and (3.2 GeV -1 ,— 4.7 GeV -1 ), following the previous analysis [20j. As given above, c 2 is 
rather well determined phenomenologically. As for C3, the value —3.4 GeV -1 is consistent 
with empirical nucleon-nucleon phase shifts, and the value —4.7 GeV -1 is the central value 
obtained from pion-nucleon scattering. There is another parameter e[(/j,), which is a com- 
bination of the 0(p 4 ) LECs and is not well known phenomenologically. Since e^(/i) is scale 
dependent, we quote its value at \x = 1 GeV in the following. 

In these formulae, the leading non-analytic quark mass dependence is given by the term 
of m\, while others (m^log(m T //x) and m%) are suppressed by additional powers of m n /Mo. 
Therefore, we also consider a simplified fit function 



M 



N 



M - 4 Cl m 2 



-ml + el(/i)m*. 



(12) 



327T/2 

When we analyse the partially quenched data set in the next section, we utilize a formula 
that is an extension of this simplified fit form. Therefore, a comparison of the simplified 
and the full fit functions (fit)]) . (jTTj) on the unitary data points provides a good test of our 
analysis. 
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I I I I I I I I I I I I I I I 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 

m n 2 [GeV 2 ] 

FIG. 2: BChPT fit of the nucleon mass for unitary points. The solid, dot, dashed, dot-dashed 
curves represent the Fit 0a, la, II, and III, respectively. 

We carry out the BChPT fits of the lattice data using these functions. The simplest 
fits are those with (fP2"|) . Since the axial-coupling g& is very well known experimentally, we 
attempt two options: (Fit 0a) a fit with fixed gA (=1.267), and (Fit 0b) a fit with gA being 
dealt as a free parameter. The fits using the 0(p 3 ) formula fflUj) are called Fit I. Again in 
this case, we attempt the fits with (Fit la) and without (Fit lb) fixing g&. For the fit using 
the 0(p^) formula ( flOi) . the lattice data do not have enough sensitivity to determine many 
parameters in the formula unless we fix gA, C2, and C3. As described above we choose gA = 
1.267, c 2 = 3.2 GeV" 2 , and c 3 = -3.4 GeV" 1 (Fit II) or c 3 = -4.7 GeV" 1 (Fit III). The 
pion decay constant is fixed at its physical value 92.4 MeV. 

We use the lattice data at five quark masses m q = 0.025, 0.035, 0.050, 0.070, 0.100. The 
data point at the smallest quark mass m q = 0.015 is not included in the fit in order to avoid 
large finite volume effect. A detailed discussion on the finite volume effect is given below. 

Figure [2] shows the ChPT fits; the resulting fit parameters are listed in Table IIHI The fit 
curves are drawn for Fit 0a, la, II and III in Figure [2] The lattice data show a significant 
curvature towards the chiral limit. All of the fit functions with fixed to the experimental 
value describe the data quite well. The results with g& a free parameter (Fit 0b and lb) 
give an important consistency check of BChPT, since the m\ term is an unique consequence 
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M [GeV] 


cifGeV" 1 ] 


e!0«) 9 A 


csIGeV" 1 ] 


cstGeV" 1 ] 


X 2 /d.o.f. 


Fit Oa 0.868(15) 


-0.97(3) 


2.89(15) [1.267] 






0.89 


Fit Ob 0.753(106) 


-1.59(56) 


6.7(3.6) 1.81(42) 






1.26 


Fit la 0.895(15) 


-0.86(3) 


3.34(16) [1.267] 






1.39 


Fit lb 0.748(104) 


-1.72(59) 


10.5(4.7) 2.13(48) 






0.88 


Fit II 0.846(13) 


-1.04(2) 


2.05(11) [1.267] 


[3.2] 


[-3.4] 


0.44 


Fit III 0.770(13) 


-1.31(2) 


1.33(12) [1.267] 


[3.2] 


[-4.7] 


0.11 



TABLE III: ChPT fit of the nucleoli mass using five unitary points m q = 0.025, 0.035, 0.050, 0.070, 
and 0.100. The values sandwiched as [• • • ] mean the input in the fit. 

of the pion loop effect in this framework. The coupling g& is in fact non-zero and roughly 
consistent with the experimental value within a large statistical error. The nucleon mass in 
the chiral limit Mq shows a significant variation, especially when the Fit III is used. 



B. Finite volume corrections 



Since the spatial extent L of the lattice is not large enough (~ 1.9 fm) for obtaining the 
baryon masses very accurately, we need to estimate the systematic error due to the finite 
volume effect. 

The finite volume correction can be calculated within BChPT, provided that the quark 
mass is small enough to apply ChPT. The nucleon mass M^{L) in a finite box of size L 3 is 



written as 



m 



M N {L) - M N {oo) = A a + A b + 0{p 5 ), 

where A a and A& represent finite volume correction at order p 3 and p 4 respectively, 

3g 2 A M ml 



A, 



A, 



m r°° I 

m7T 2 f fJ o dx^K (L\n\^/M^ + ml(l-x)) 
— 2_^ (2ci-c 3 ) — — + c 2 - 



3m 



K x {L\n 




L 


n 





(L|n|ra,r) 



(13) 

(14) 
(15) 



Here, the functions ^0(^)5 K\(x) and K 2 {x) are the modified Bessel functions, which asymp- 
totically behave as exp(— x) for large x. The sum runs over a three dimensional vector n of 
integer components, and |n| denotes a/^ 2 ". 

In the following we make the following two different analyses for the finite volume effect. 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 
m n 2 [GeV 2 ] 

FIG. 3: Chiral fit of the corrected data (diamonds). Solid and dashed curves represent the fits using 
5 and 6 heaviest data points, respectively. For a reference, we also show the raw data (circles). 

M [GeV] dtGeV 1 ] gQz) X 2 /d.o.f. 

Fit 0a(5pt) 0.793(15) -1.04(3) 2.68(15) 1.86 
Fit 0a(6pt) 0.808(13) -1.02(3) 2.79(13) 1.82 

TABLE IV: Chiral fit parameters for the finite volume corrected lattice data. Results using all 6 
data points "Fit 0a (6pt)" and those using 5 heaviest data points "Fit 0a (5pt)" are listed. The 
fit function is (|12p with a fixed axial coupling qa = 1.267. 

1. We correct the data for the finite volume effect using the above formula. For the input 
parameters M , Qa, c% (i = 1-3), we use the nominal values (M = 0.87 GeV, qa = 
1.267, ci = -1.0 GeV" 1 , c 2 = 3.2 GeV~\ and c 3 = -3.4 GeV" 1 ). The size of the 
finite volume corrections varies from —0.3% (heaviest) to —4.0% (second lightest) and 
—6.5% (lightest). The chiral fit is then made for the corrected data points using the 
simplified 0(p 3 ) formula ffT2l for heaviest 5 or 6 heaviest data points. The result is 
shown in Figure [3] and the fit parameters are listed in Table IIVI After correcting the 
finite volume effect, there is a 5-8% decrease in Mq and 4-7% increase in the magnitude 
of the slope in the chiral limit |ci|. The results of the fits with 5 or 6 data points are 
consistent with each other. 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 
m n 2 [GeV 2 ] 



FIG. 4: Chiral fit with the 0(p 4: ) formula plus the finite volume effect A a + Dashed, solid and 
dotted curves represent the fit results in the infinite volume from the Fit la, II and III, respectively. 
For a reference, we show the finite volume corrected data points (diamonds). 

2. We fit the data with the fit functions including the finite volume corrections, i.e. at 
0(p 3 ) the function is flED plus A a (Fit II); at C(p 4 ) the function is (El) plus A a + A b 
(Fit III) . Figure H] shows the fit curves after subtracting the finite volume piece A a or 
A a + Ab, which consistently run through the finite volume corrected data points. The 
fit parameters are listed in Table |Vj We find that after taking the finite volume effect 
into account M decreases by 3-8% and |ci| increases by 9%. The 5-point and 6-point 
fits are consistent with each other within two standard deviations. 

Comparing the fit parameters obtained with (Tables IIVI and |V]) and without (Table IIIB the 
finite volume corrections, we observe that the deviation due to the finite volume effect is 
smaller than the uncertainty of the fit forms. 

There is also a finite volume effect due to fixing the topological charge in our simulation. 



This can be estimated using ChPT as in 



33 



34J. The estimated corrections are fairly 



small, —(0.3-0.7)% depending on the quark mass. Compared to the statistical error and the 
conventional finite volume effect, we can safely neglect the fixed topology effect. 



15 



M 




9a 


C2 


C3 


X 2 /d.oi. 


[GeV] 


[GeV" 1 ] [GeV- 3 ] 




[GeV" 1 ] 


[GeV" 1 ] 




Fit Ia(5pt) 0.852(15) 


-0.90(3) 3.18(16) 


[1.267] 






1.81 


Fit Ia(6pt) 0.870(13) 


-0.88(2) 3.31(13) 


[1.267] 






1.86 


Fit II(5pt) 0.778(12) 


-1.08(2) 1.70(13) 


[1.267] 


[3.2] 


[-3.4] 


1.19 


Fit II(6pt) 0.794(10) 


-1.06(1) 1.83(11) 


[1.267] 


[3.2] 


[-3.4] 


1.76 


Fit III(5pt) 0.694(12) 


-1.35(2) 0.84(14) 


[1.267] 


[3.2] 


[-4.7] 


0.24 


Fit III(6pt) 0.723(10) 


-1.32(1) 1.10(12) 


[1.267] 


[3.2] 


[-4.7] 


3.37 



TABLE V: Results from the chiral fit including the finite volume corrections. The finite volume 
effects are included to A a for the Fit la (at 0(p 3 )) and to A a + A b for the Fit II and III (at C(p 4 )). 





w/o FVCs 


w/ FVCs 




5 pt 


5 pt 6 pt 


Fit 0a 


52.2(1.8) 


56.7(1.8) 55.1(1.5) 


Fit la 


45.1(1.7) 


48.9(1.7) 47.2(1.5) 


Fit II 


56.5(1.2) 


59.5(1.2) 58.2(1.0) 


Fit III 


71.8(1.2) 


75.1(1.2) 72.7(1.0) 



TABLE VI: Nucleon sigma term ovjv [MeV] with and without the finite volume corrections (FVCs). 
C. Nucleon sigma term 

Using the fits in the previous subsections we obtain the nucleon sigma term by differen- 
tiating the nucleon mass with respect to the quark mass as 

dM N 



E 



q=u,: 



(16) 



m q =m ud 



Since the value of the physical up and down quark mass is very small, we may extract the 
physical value using the leading-order ChPT relation 

o dM N 



dml 



(17) 

771^=135 MeV 

Table [VI] shows the results from the several fit forms with and without the finite volume 
corrections (FVCs). 
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Due to the curvature observed in Figures EHU that is largely explained by the non-analytic 
term in the BChPT formulae, g^n is enhanced toward the chiral limit. Compared with 
the value at around the strange quark mass, a n N is about three times larger, depending on 
the details of the fit ansatz. 

The largest uncertainty comes from the chiral extrapolation. In fact, the Fit III gives 
significantly larger value of a n N than those of other fit ansatz. It is expected from the plot 
of chiral extrapolation, Figure [2], where the Fit III (dot-dashed curve) shows a steeper slope 
near the chiral limit. The finite volume effect is a sub-leading effect, which is about 9%. We 
take the Fit Oa (gA fixed, FVCs not included) as our best fit, and take the variation with fit 
ansatz and FVCs as an estimate of the systematic error. We obtain 

a nN = 52(2) stat (l 2 °) extrap (^) FVE MeV, (18) 

where the errors are the statistical and the systematic due to the chiral extrapolation (extrap) 
and finite volume effect (FVE). This result is in good agreement with the phenomenological 
analysis based on the experimental data at the Cheng-Dashen point a^N — 55 ~ 75 MeV, 
which is discussed in the Introduction. 

V. ANALYSIS OF THE PARTIALLY QUENCHED DATA POINTS 
A. Fits with partially quenched ChPT formula 

As described in Section [Til the partial derivatives in terms of the valence and sea quark 
masses, m va i and m sea respectively, are necessary in order to extract the connected and 
disconnected-diagram contributions separately, hence to obtain the strange quark content 
y defined in (J2J). It is possible with the lattice data in the so-called partially quenched 
set-up, i.e. the valence quark mass is taken differently from the sea quark mass. Since the 
enhancement of ovjv towards the chiral limit is essential for reliable determination of the 
nucleon sigma term, we should use the chiral perturbation theory formula for baryons in 
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partially quenched QCD, which is available for two-flavor QCD [36l. 137] . At 0(p 3 ), it reads 



M N = B m + B l0 {m v :f + B m {mi s f 



167T/2 I 12 

-f| [-19(mD 3 + 10(m-) 3 + 9mr(mr) 2 )] 

[-13(mD 3 + 4(r<) 3 + 9m™ (r<) 2 )] } 



91 
12 



3 

+B 20 (mm 4 + 5nW) 2 Kf + 5o 2 (mr) 4 , (19) 

where m™, m^ s , and m^ s denote the pion mass made of valence- valence, valence-sea, and sea- 
sea quark combinations, respectively. At this order of the chiral expansion, one can rewrite 



this formula in terms of m va i and m sea using the leading-order relations (m™) = 2Bm 



val) 



vs \ 2 



-B(m va i + m sea ), and {m s J) 2 = 2Bm seau . The parameter B is determined as B 



= 1.679(4) GeV through the ChPT analysis of pion mass 35|. The coupling constant 
gA represents the nucleon axial charge as before, and g\ is another axial-vector coupling 
characterizing the coupling to the rj meson. They are related to the standard F and D 
parameters of BChPT as gA — F + D and g\ = 2(F — D). Numerically, the values of F 
and D are obtained from the hyperon decay as F = 0.52(4) and D = 0.85(6) (see [38], for 
instance), which imply g\ = —0.66(14). In the following, whenever we need nominal values 
of gA and gi, we set = 1.267 and g\ = —0.66. 

Strictly speaking, there are also contributions from the decuplet baryons. In our analysis 
we have integrated out the Delta resonance and expand the contribution in terms of (m^/A) 2 
with A = ira — Mjy. Then these contributions can be absorbed into the analytic terms in 
(EH}. 

We fit the quark mass dependence of the nucleon mass with the partially quenched ChPT 
formula (fT9l . The independent fit parameters are B 00 , £> i, B w , B n , B 2 o, B 02 , <7i, and gA- 
Instead of making all these parameters free, we also attempt a fit with fixed gA and g\ (Fit 
PQ-a), a fit with fixed gA (Fit PQ-b). The fit with all the free parameters is called the Fit 
PQ-c. 

Figure [5] demonstrates the result of the partially quenched ChPT fit. It shows the sea 
quark mass dependence at eight fixed valence quark masses. Data are nicely fitted with the 
formula ([191 . The fit results are listed in Table IVHl All the parameters are well determined 
except for the term B^m^) 4 , for which the data do not have enough sensitivity. Finite 
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-^00 Bqi B w Bu B 2 o B 02 

[GeV] [GeV" 1 ] [GeV" 1 ] [GeV~ 3 ] [GeV~ 3 ] [GeV~ 3 ] 



9l 



g A x 2 /d.o.f. 



Fit PQ-a 0.87(2) 0.47(10) 3.37(3) -0.94(2) 3.77(2) 0.17(15) [-0.66] [1.267] 1.82 
Fit PQ-b 0.86(2) 1.13(11) 2.71(4) 0.97(8) 1.81(11) 0.25(15) -0.378(14) [1.267] 1.28 
Fit PQ-c 0.92(4) 0.76(23) 1.98(39) 0.43(31) 0.95(43) -0.03(23) -0.29(5) 0.93(22) 1.28 

TABLE VII: Fit results with the partially quenched ChPT formula. The values sandwiched as 
[• • • ] mean the input in the fit. 




• 


am 

V 


=0.025 

ii 


■ 


am 

V 


=0.035 

ill 


♦ 


am 

V 


=0.050 

ill 


A 


am 

V 


=0.060 

ill 


4 


am 


^=0.070 




am 


^=0.080 




am 


=0.090 

ill 


* 


am 


^=0.100 



0.3 0.4 0.5 

(m "f [GeV 2 ] 



0.6 



0.7 



0.8 



FIG. 5: Partially quenched nucleon masses and fit curves (Fit PQ-b). 

volume corrections are not taken into account. 

By reducing the parameters to the case of the unitary point m val = m sea (m = -B005 
ci = —(Bqi + Bio)/ 4, ei = B20 + Bn + -B02), it is easy to see that the results from Fits 
PQ-a and PQ-b are consistent with the Fit 0a for the unitary points. Figure [6] shows the 
reduction to the unitary point. The value of the nucleon sigma term u^n obtained from this 
reduced fit parameters is 53.3(1.8), 53.2(1.9) and 41.3(6.6) MeV for the Fits PQ-a, PQ-b, 
and PQ-c, respectively. These values are in good agreement with our analysis of the unitary 
points (US]). 

Another important observation from Table IVHI is that the Fit PQ-c, for which is a 
free parameter, gives much better constrained g& than the Fit 0b of the unitary points. 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 
m n 2 [GeV 2 ] 

FIG. 6: Result of the PQChPT fit for unitary points. Solid, dashed and dotted curves represent 
the fit results from the Fit PQ-a, PQ-b and PQ-c, respectively. 

This is because the partially quenched analysis uses much more data points: 40 data points 
compared to 5 in the unitary case. It is remarkable that both and g\ can be determined 
with reasonable accuracy. 

B. Sea quark content of the nucleon 

Once the valence and sea quark mass dependence is identified using the formula (fl~9j) . we 
can obtain the partial derivatives with respect to m va \ and m sea to obtain the connected and 
the disconnected contribution to the nucleon sigma term cr nN as defined in © and (jSJ). 

Figure [7] shows the partial derivatives with respect to m va i (left panel) and to m sea (right 
panel) evaluated at the unitary points m va i = m sea . In the plots, the fit results are plotted as 
a function of m va i = m sea . For both contributions, we clearly find an enhancement towards 
the chiral limit. Results with different fit ansatz show slight disagreement near the chiral 
limit, which indicate the size of the systematic uncertainty. 

Numerical results at the average up and down quark mass and at the physical strange 
quark mass are listed in Table IVIIIl The values in the lattice unit am u d = 0.0034(1) and 
am s = 0.084(2) are determined from a partially quenched analysis of the meson spectrum 
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am . = am 

val sea 



FIG. 7: Connected (left) and disconnected (right) contributions to the nucleon sigma term eval- 
uated at m va i = m sea . Solid, dashed and dotted curves represent the results from the Fit PQ-a, 
PQ-b and PQ-c, respectively (thick lines). The error curves are represented by the thin lines. 





dM N dM N 
am i a a 


dM N dM N 


Fit PQ-a 
Fit PQ-b 
Fit PQ-c 


0.0034 7.92(8) 1.47(32) 
0.0034 6.68(8) 2.72(33) 
0.0034 5.27(75) 1.99(50) 


0.084 2.75(3) 0.28(14) 
0.084 2.84(3) 0.28(14) 
0.084 2.81(3) 0.26(14) 



TABLE VIII: Connected and disconnected contributions to the nucleon sigma term, evaluated 
at the average up and down quark mass am q = 0.0034 and at the physical strange quark mass 
am q = 0.084. 



35|. 

Figure M shows the ratio of the disconnected and connected contribution to the sigma 
term (N\(uu + dd)\N) disc / (N\(uu + dd)\N) conn evaluated at the unitary points m sea = m va j. 
We find that the sea quark content of the nucleon is less than 0.4 for the entire quark mass 
region in our study, so that the valence quark content is the dominant contribution to the 
sigma term. This is in striking contrast to the previous lattice results in which the sea quark 
content equal to or even larger than the valence quark content was found. 
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FIG. 8: Ratio of the disconnected and connected contribution to the sigma term for unitary points 
(m sea = ^vai)- Solid, dashed and dotted curves represent the results from the Fit PQ-a, PQ-b and 
PQ-c, respectively (thick lines). The error curves are represented by the thin lines. 

C. A semi-quenched estimate of the strange quark content 



Rigorously speaking, it is not possible to extract the strange quark content (iV|ss|iV) 
within two-flavor QCD. The problem is not just the strange quark loop is missing, but it 
is not possible to evaluate the disconnected contribution at the strange quark mass while 
sending the sea and valence quark masses to the physical up and down quark mass with the 
partial derivatives within (partially quenched) two-flavor QCD. For the final result, therefore, 
we should wait for a 2+1-flavor QCD simulation, which is in progress 40J. Instead, in this 
work, we provide a "semi-quenched" estimate of the strange quark content assuming that the 
disconnected contribution gives a good estimate of the strange quark effect when evaluated 
at the strange quark mass for both m val and m sea . 

We define our semi-quenched estimate of the parameter y as the ratio of the strange 
quark content (disconnected contribution at m va \ = m sea = m s ) to the up and down quark 
contributions (connected plus disconnected contributions at m va i = m sea = m ud ) following 
Ref. [la ] . Taking the result from the Fit PQ-b as a best estimate, we obtain the parameter 
V as 

y Nf=2 = 0.030(16) stat (^) cxtrap (^) ms , (20) 
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where the errors are statistical, the systematic errors from chiral extrapolation and from the 
uncertainty of m s . The chiral extrapolation error for the strange quark content is estimated 
by the differences of the results of Fit PQ-a,b and c, while that for the up and down quark 
content is estimated by the differences of the results of Fit 0, I, II and III. We also note that 
there may be an additional ~10% error from finite volume effect as discussed in Section HVl 
but it is much smaller than the statistical error in our calculation. 



VI. DISCUSSION 

We found that the disconnected contribution to the sigma term is much smaller than 
the previous lattice calculations with the Wilson-type fermions y ~ 0.36 ~ 0.66 [l3, [l4, Ijj] 
(except for 16j as explained below). The authors of 16| found that the naive calculation with 
the Wilson-type fermions may over-estimate the sea quark contents due to the additive mass 
shift and the sea quark mass dependence of the lattice spacing. The key observation is that 
the additive mass shift is large depending significantly on the sea quark mass. Therefore, in 
order to obtain the derivative ([6]) one must subtract the unphysical contribution from the 
additive mass shift. This problem remains implicitly in the quenched calculations, since the 
derivative must be evaluated at the value of the valence quark mass even when the sea quark 
mass is sent to infinity. (There is of course the more fundamental problem in the quenched 
calculations due to the missing sea quark effects.) 

Another problem is in the conventional scheme of setting the lattice scale in unquenched 
simulations. In many dynamical fermion simulations, the lattice spacing is set (typically 
using the Sommer scale ro) at each sea quark mass, or in some cases, the bare lattice coupling 
(3 is tuned to yield a given value of r independent of the sea quark mass. This procedure 
defines a renormalization scheme that is mass dependent, because the quantity ro could have 
physical sea quark mass dependence. Since the partial derivative (jEJ) is defined in a mass 
independent scheme, i.e. the coupling constant does not depend on the sea quark mass, one 
has to correct for the artificial sea quark mass dependence through r when one calculates 
the nucleon sigma term. Combining these two effects, the authors of [if]] found that their 
unsubtracted result y = 0.53(12) is substantially reduced and becomes consistent with zero: 
y = —0.28(33). The conclusion of this analysis is that the previous lattice calculations giving 
the large values of y suffered from the large systematic effect, hence should not be taken at 
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their face values. 

Our calculation using the overlap fermion is free from these artifacts. The additive mass 
shift is absent because of the exact chiral symmetry of the overlap fermion. The lattice 
spacing is kept fixed in our analysis at a fixed bare lattice coupling constant. We confirmed 
that this choice gives a constant value of the renormalized coupling constant in the (mass 
independent) MS scheme through an analysis of current-current correlators [39]. Therefore, 
the small value of y obtained in our analysis fl20|) provides a much more reliable estimate 
than the previous lattice calculations. 



VII. SUMMARY 



We study the nucleon sigma term in two-flavor QCD simulation on the lattice with exact 
chiral symmetry. Fitting the quark mass dependence of the nucleon mass using the formulae 
from Baryon Chiral Perturbation Theory (BChPT), we obtain a n N = 53(2) (t 2 )) MeV, 
where our estimates of systematic errors are added in quadrature. This is consistent with 
the canonical value in the previous phenomenological analysis. Owing to the exact chiral 
symmetry, our lattice calculation is free from the large lattice artifacts coming from the 
additive mass shift present in the Wilson-type fermion formulations. 

We also estimate the strange quark content of nucleon. From an analysis of partially 
quenched lattice data, we find that the sea quark content of the nucleon is less than 0.4 for 
the entire quark mass region in our study. The valence quark content is in fact the dominant 
contribution to the sigma term. Taking account of the enhancement of (N\(uu + dd)\N) 
near the chiral limit, the parameter y is most likely less than 0.05 in contrast to the previous 
lattice calculations. 

By directly calculating the disconnected diagram we may obtain further information. For 
instance, the effect of the strange quark loop on the dynamical configurations with light up 
and down quarks can be extracted. Such a calculation is in progress using the all-to-all 
quark propagators on the lattice. Another obvious extension of this work is the calculation 
including the strange quark loop in the vacuum. Simulations with two light and one strange 
dynamical overlap quarks are on-going 40]. 
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